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Abstract 

In this paper, a study on discrete-time coined quantum walks on the 
line is presented. Clear mathematical foundations are still lacking for 
this quantum walk model. As a step towards this objective, the following 
question is being addressed: Given a graph, what is the probability that 
a quantum walk arrives at a given vertex after some number of steps? 
This is a very natural question, and for random walks it can be answered 
by several different combinatorial arguments. For quantum walks this 
is a highly non-trivial task. Furthermore, this was only achieved before 
for one specific coin operator (Hadamard operator) for walks on the line. 
Even considering only walks on lines, generalizing these computations to 
a general 5*17(2) coin operator is a complex task. The main contribution 
is a closed-form formula for the amplitudes of the state of the walk (which 
includes the question above) for a general symmetric SU{2) operator for 
walks on the line. To this end, a coin operator with parameters that alters 
the phase of the state of the walk is defined. Then, closed-form solutions 
are computed by means of Fourier analysis and asymptotic approximation 
methods. We also present some basic properties of the walk which can 
be deducted using weak convergence theorems for quantum walks. In 
particular, the support of the induced probability distribution of the walk 
is calculated. Then, it is shown how changing the parameters in the coin 
operator affects the resulting probability distribution. 
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1 Introduction 



The design of quantum algorithms is nowadays one of the major problems in the 
quantum computing community. The strategies for writing classical algorithms 
as divide and conquer, dynamic programming, etc, are not easily adapted to the 
quantum paradigm. Strategies for designing quantum algorithms are phase am- 
plification, phase estimation, to name a few. As an example of the applications 
of these strategies, Grover's algorithm uses the amplitude amplification tech- 
nique, and Shor's algorithm relies in reductions to order finding and phase esti- 
mation [T7|. Therefore, it becomes necessary the study of different approaches 
to improve the efficiency of the search. 

One of the emergent alternatives for the design of algorithms are quantum 
walks, in direct analogy to random walks in classical computing. Random walks 
showed to be a successful tool for designing algorithms, and the same success 
is expected in the quantum paradigm. Results in this field showed that quan- 
tum walks can outperform its classical counterpart by exploiting quantum me- 
chanical effects such as interference and superposition, giving and exponential 
speedup for certain types of graphs, and polynomial speedup for some practical 
applications [Tl [TO]. 

There are two types of discrete-time quantum walks. Quantum Markov 
Chains and Coined Quantum Walks. This paper is about the latter, and from 
now on we will refer to this model simply as quantum walk when is obvious 
from the context. 

The field of quantum walks is very recent, and still lacks a solid mathemat- 
ical foundation. Markov chain quantum walks already started to build these 
foundations by establishing a direct connection to classical Markov chains using 
algebraic techniques [IS]. However, coined quantum walks are not having the 
same luck, and it seems that mathematical techniques for random walks simply 
do not work. 

1.1 Related Work 

Coined Quantum Walks are defined by the application of two unitary operators 
S and C, where C (coin operator) decides which vertex to move onto, and S 
(shift operator) performs the actual movement of the walk given the direction 
decided by C. Ambainis [1], Kempe [10] and Konno [H] give good surveys 
of this model. There are several studies of this walk for specific graphs. On 
the line, Ambainis et al. [3] and Chandrashekar, Srikanth, and Laflamme [6] 
show that the variance of the induced probability distribution has a quadratic 
improvement over the classical walk (i.e. for t steps, V = 0{t^) and classically 
V = 0{t)). Konno computed the induced probability distribution using path 
integrals [H] and via a weak limit theorem [T3]. In the hypercube, Kempe [TT] 
shows that the hitting time from one corner to its opposite is exponentially 
faster, while Moore and Russell [15] gives the same speed-up for the mixing 
time. For practical applications there are algorithms for hypercubes and grids. 
For the hypercube, Shenvi et al. [20] gives an algorithm for solving SAT with a 
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quadratic improvement, while Potocek et al. [TH] gives an improvement of the 
same algorithm on the success probability. For grids, Ambainis et al. [4] show 
a quadratic speed-up and presents a general framework for analyzing quantum 
walks. Also, Ambainis [5] gives an optimal algorithm for element distinctness 
over the Jhonson graph with a quadratic speed-up. 

Quantum walks on the line is probably the most studied quantum walk 
model. Interest on this matter started in computer science with Ambainis, 
Bach, Nayak, Vishwanath, and Watrous [3] , where notions of hitting and mixing 
times were introduced. In the same piece of work, they computed a closed-form 
formula for the induced probability distribution of a Hadamard walk (i.e. a 
quantum walk with a Hadamard operator as coin). Furthermore, their formula 
gives a complete characterization of the amplitudes in the state of the walk in 
the asymptotic limit. 

It is known that the dynamics of the walk is controlled by the coin operator 
[TU] . Thus, depending on the application, a good choice of the coin could make a 
great difference. This motivated the study of quantum walks on the line moved 
by a general SU (2) operator, which has four independent variables. However if 
we consider only the resulting probability distribution, one variable is enough; 
i.e. any probability distribution resulting from a quantum walk on the line can 
be simulated by a general rotation around the z axis with parameter 9. Nayak 
and Vishwanath (TB] gave an intuitive description of the probability distribu- 
tion based on the stationary phase method without giving an explicit formula 
for it, and without considering the amplitudes of the state of the walk. Chan- 
drashekar, Srikanth, and Laflamme studied generalized walks using a SU (2) 
coin operation. They present an approximate formula for the amplitudes of the 
state of the walk. However, their results were based in numerical experiments 
rather than a complete analytically deducted formula. Grimmet, Janson, and 
Scudo [8] showed a ballistic spreading of the walk and they gave an expression 
for the limit distribution using weak convergence theorems. 

1.2 Contributions 

As a step toward finding mathematical foundations of quantum walks, in this 
paper the following question is being addressed: Given a graph, what is the 
probability that a quantum walk arrives at a given vertex after some number 
steps? This is a very natural question, and for random walks it can be answered 
by several different combinatorial arguments [22]. 

The main contribution of this paper is a closed- form formula^ for the question 
above for a general symmetric SU{2) operator for walks on the line (Theorem 
[3]). Furthermore, the formula characterizes the amplitudes of the state of the 
walk in the asymptotic limit. In comparison to the previous works mentioned 
before (Nayak and Vishwanath jl6], Chandrashekar et al. [6]), the closed- form 
formulas derived in this paper were analytically computed for the amplitudes 

quantity /(n) is in closed-form if we can compute it using at most a fixed number of 
"well-known" standard operations, independent of n [7]. 
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Tabic 1: Known results for different coins for walks on the line. 



Coin 


Amplitudes of the 
state 


Probability distribu- 
tion 


Hadamard 


closed-form [16^ 


closed-form [H [H] 


SU{2) 


numerical results [6] 


numerical results [6] , 
closed-form [TQ 


Symmetric 
SU{2) 


closed-form [this work] 


closed-form [this work] 


U{2) 


explicit formula (not 
closed-form) [12 [13] 


explicit formula (not 
closed-form) [l3l 



of the state of the walk (including the induced probability distribution) for a 
symmetric SU{2) operator (Table [1] shows more clearly these differences). Also, 
in a seminal work, Konno \12\ I13j gave explicit expressions for the amplitudes 
of a U[2) coin, using a discrete path integral method in a clever way. However, 
these expressions were not in closed-form, as we claim in this work. Further- 
more, we show how to compute the errors in the asymptotic approximation, 
something that was missing from previous works in the literature. To this end, 
in Section [2] a coin operator with parameters that alters the phase of the state of 
the walk on the line is proposed. The coin operator is inspired by the quantum 
algorithm for SAT proposed by Hogg [H]. In that work, in order to implement 
heuristics for quantum algorithms, the author proposed to add parameters to 
the unitary operation of a search algorithm. This way, the situation is similar 
to classical algorithms where a tunable set of parameters are adjusted according 
to the problem. After defining the coin operation, we compute the spectrum 
of the unitary evolution operator of the walk using Fourier analysis. In Section 
[3J after having obtained the eigenspectrum of the walk, we apply the inverse 
Fourier transform to obtain the state of the walk in terms of Fourier coefficients. 
To compute a closed-form solution in the asymptotic limit from the Fourier co- 
efficients, we applied the Euler-Maclaurin formula [5] and the steepest descent 
method for asymptotic approximation of integrals |23| . This method is in fact 
stronger than the stationary point method from |16j and [21] . where the au- 
thors use it to study the asymptotics of the resulting probability distribution 
from coin operators with real eigenvalues. With the steepest descent method we 
can compute the amplitudes of the state of the walk resulting from any complex 
unitary operator. In Section 21 we compute the error terms for the approxima- 
tions made, which can be derived from the employed methods. Finally, some 
basic properties of the walk are examined by means of weak convergence the- 
orems [H]. The support of the induced probability distribution of the walk is 
computed, and then we argue how changing the parameters in the coin operator 
affects the resulting probability distribution. 
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2 Quantum Walks with Phase Parameters 

In this section, a coin operation with parameters is proposed. Then, using 
Fourier analysis, integral-forms for the amplitudes of the walk on the line are 
computed. Later in Section[31 it is shown how to solve these integrals and derive 
a closed-form in the asymptotic limit. 

2.1 Walks on the Line with Phase Parameters 

Here we define quantum walks on the line, and introduce the coin operator used 
in this research. 

Definition 1. Let He = span{|-s— ) ,|— >)} and Hs = span{\n) : n e Z}. The 
state of the walk Y^t) = X^n IV^tl*^)) '^^ time t is defined over the joint space 
He <8) Hs with basis states {\d,n) : \d) € Hc,\n) G Hs}, where liptin)) = 
a^(n)|d, n) and af{n) is the amplitude at time t in direction d and posi- 
tion n. ^^.so „ = I. 

For the analysis of the walk on the line we consider the projection at time t 
onto position n as a 2 dimensional vector, i.e. 



with (n) and (n) representing the amplitude of the walker at position n 
at time t going left and right respectively. The probability of being at position 
n at time t is thus given by 



Throughout the paper, the initial condition is considered as \tpo{0)) ~ [acTi Q^cT]"^ 
and IV'oN) = [0,0p for n ^ 0, with \a^\^ + \a^\^ = 1. 

The quantum walk is defined by the way it moves at each time step. This is 
captured by the following definition. 

Definition 2. The time evolution of the walk is given by 

|*f) = U l^t^i) , or equivalently, |*f) = [/* |*o> , 

where U — S{C ® I) is a unitary operator defined on the Hilbert space of the 
whole system He ® Ha, I is the identity matrix acting on Hs, C is the coin 
operator acting solely on He, and S is the shift operator in charge of performing 
the walk. 

According to this definition, the walk first choses a direction of movement 
using C, and then moves with operator 5*. In order to move, operator S needs 
to be conditioned on the coin space in the following way. 



P,{n) ^ \{Mn)\Mn))\' = latin)]' + latin)]'. 



(1) 





n 



5 



flood) 




Figure 1: Quantum walk on the line with different values of phase parameters. 
The variance of the walk changes depending on ti and T2 . Since the probabilities 
at odd positions are 0, those points are not plotted. 



Definition 3. The coin operator is defined by C = HTH, where H is the 
Hadamard operato^ in charge of mixing amplitudes among states, and T = 
gjiTTii ^)(^ I + e""^^! —!>)(—> I is the diagonal phase adjustments with ti,T2 € 
[0,1]. 

Let a = e'^'^i + e*'^'^^ and b = e^'^'^i — 6"^"^^ ^ xhcn, the resulting operator can 
be written as 

which have the following effect on He 

|H^(l/2)a| H + (1/2)61 
1^)^(1/2)61 H + (l/2)aH). 

Figure [1] shows the dynamics of a walk using C as coin. For different val- 
ues of the phase parameters ti and T2 the variance of the induced probability 
distribution changes. 

The state of the walk at time t can be related to the state at time t + 1 
according to the following lemma. 

Lemma 1. 

\^t+i{n)) = M+\Mn - 1)) + AUMn + 1)) (3) 
^The Hadamard operator is defined as H = -y= 



a b 
b a 
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whe 



M. = 



■ 

(1/2)6 (l/2)a 



id A/_ = 



■(l/2)a (1/2)6- 




Proof. Let l^ff) = a^(n)| , n) +al^{n)\ — >, n) be the state at time t. Also 
denote the amphtudcs after applying operators C and S as 

(C® /) I*,) = ^«r(n)'| +ar(n)'| ^,n), 



Now let 1^1/4+1 



be the state at time t + 1. The amplitudes of this 



state are related to the amplitudes of j'l't) in the following way 







■<(")"■ 




ar(n + l)'' 













The contributions to the amplitudes of state |5't+i) come from position n + 1 
for the upper component, and from n — 1 for the lower component by definition 
of operator S. The amplitudes corresponding to the state after applying C are 
computed as follows: 



- _{l/2)hat{n 



-l) + {l/2)ha^{n + l) 
l) + (l/2)aar(n+l) 



and the same for C\ijjt{n— 1)). Thus 
\'4^t+i{n)) 



(l/2)aat{n + 1) + (l/2)6at-^(?i + 1) 
(l/2)6ar(n-l) + (l/2)aar("-l) 
= M+\Mn-l))+M4i^t{n + l)), 



where 





(1/2)6 (l/2)a 



and Af_ 



■(l/2)a (1/2)6 




□ 



2.2 Analysis 

One approach to the analysis of quantum processes is the path integral approach. 
This method explicitly computes the amplitude of a certain state as the sum 
over all possible paths leading to that state [3l [12] . Solving a path integral is 
known to be hard, and we avoid this by following the steps of [31 [TSl [TT] known 
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as the Schrddinger approach. Given the translational invariance of the walk, it 
has a simple description in Fourier space [3] . The Fourier transform of the walk 
is analyzed and then transformed back to the original domain. 

The quantum Fourier transform |17l of a wave equation is defined by 



and the corresponding inverse Fourier transform is then 



\Mn)) = I 



— ikvi 



Mk))dk 



(4) 



(5) 



Applying gl) to JS]) we get 

71 

= e'^M+ e''^-("-i) |V;t(n - 1)) 

n 

+ e-*'=A/_ Y 6'''^"+^' I Vt (ri + 1)) 



Mk) 



Aik) 



= e"'M+ 

Then, the time-evolution in Fourier space is given by 

i^t+i{k))=Mk Mk) 



where = e^^Mj^ + e ^^M-. In matrix form 



be 



■,k ^^ik 



(6) 



(7) 



In general, the state at time t is given by the i-th power of operator Mk 
applied to the initial state 

Mk)) = Mi \Mk)) ■ (8) 
The following lemma shows the eigenspectrum of operator Mk ■ 

Lemma 2. Let Mk he a unitary matrix as in The eigenvalues and eigen- 
vectors of Mk are 

Xj {k) ^1/2 {a cos k ± ^/iP -d^ sin^ fc) 



and 



\\,{k)) = N,{k) 



-iasinfci \/\P- — sin^ k 
be"^ 
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respectively, with j = 1,2. Furthermore, Nj{k) is a normalization coefficient 
given by 



N,{k) = 



—ai sin k ± \/ h"^ — a? sin^ k 



-1/2 



Proof. The characteristic polynomial of Mk is determined by det{Mk — A/) = 0. 
Then 

^2 ^2 

det{Mk - XI) = - aAcosfc + — - — . 

Solving the equation gives the eigenvalues 

a cos k ± \/\P- — ci^ sin^ fc 
Aj(fc) = ^ . 

for J = 1, 2. In order to find the eigenvectors, we solve the following system of 
linear equations 



(Mfe-A,(fc)/) 











et a 





de-'*^ -Aj-(fc)) +2/^-^6-*'= 



b ik 



X,{k)) 





0" 








this vector is still an eigenvector, multiply yj and by be and obtain 



Xj 




■-ae*''' + 2A/ 




—ai sin k ± \/b'^ — a'^ sin^ fc 













Then Nj{k) is 1 divided by the f2-norm of this vector, and multiply the eigen- 



vectors by Njik) to normalize them. 
Diagonalize ([7]) to obtain 



□ 



Ml^ J2 A,(fc)*|A,(fc))(A,(fc)|, 
je{i,2} 

where Ai(fc) and A2(fc) are the eigenvalues with corresponding eigenvectors 
|Ai(A:)) and |A2(fc)). Now apply the diagonalized operator to the time evolu- 
tion ([HI and obtain the following form 

j 

= J2{x,{k) Mk)) X,{kY \X,{k)) . (9) 



The initial state is [a^,a(/'] , 



and in Fourier space becomes \'>po{k)) 



[a^,a^]'^ for all k E [— 7r,7r]. To write equation ([9]) in a simpler way, define 
0(fc) = (A,(fc)|^o(fc)) 

= a^Njik) (-msinfci A/52^r^2~^y _^Q,^^^.(fc)^*g-*fc^ (10) 
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where * is the complex conjugate. This can be expressed in matrix form as 



"el(fc)' 









{-ia sin k+Vb"^ - a? sin^ k)* h*e-'^Ni{k) 
{~ia sin k- \/b^ - a? sir? k)* h*e-'^N2{k) 

The state of the walk at time t can be expressed by 



Ml 



Mk)) = Y.\]{k)^j{k)\\,{k)) . 



(11) 



Let a^(fc) and aj'ik) be the amplitudes of the state i^t{k)j in Fourier space 

going left and right respectively. Then, by equation PT|) and Lemma [5] these 
amplitudes are 

at{k) = ^\j{ky^j{k)Nj{k) (-msinfczb -a2 sin^ fcj (12) 



and 



a^{k)=Y.^,{kr^,{k)N,{k)be'\ 



(13) 



The final step is to reverse back to the original domain of the walk. This is 
done by applying ^ to and (fT3|) . 



-!- / ^S,,{k)Nj{k)\]e-'^'' (-msinfc ± y/b^ - a? sin^ k) dk (14) 



and 



«r W = ^ fY.^Ak)NAk)be''\]e-^''"dk. (15) 



Note that a discrete walk is being approximated by an integral. The Euler- 
Maclaurin summation formulgfl gives the error term for these approximations. 

Equations and (flSl) can be solved by the steepest descent method from 
complex analysis, obtaining this way closed-form solutions. This is done in the 
next section. 



3 Asymptotic Approximation 

In this section it is shown how to find close-form solutions to the integrals (|14p 
and ([TS]) . First, in Section [3.11 the technique used in this research known as 
the steepest descent method is briefiy explained. Then, in Section [X^ the same 
technique is applied to the integral- forms of the walk (equations and ([T5)) ). 

'EL. /(") = /a fi^)d^ + i^^^^ + Er=i (i^ (b) - /(^''-D (a)), where eaeh 

Bi is a Bernoulli number [5]. 
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3.1 Steepest Descent Method 

Here one of the most powerful methods for asymptotic approximation of inte- 
grals is briefly explained. The method is known as Steepest Descent Approxi- 
mation or Saddle Point Method. For a deeper understanding on this technique 
refer to 

The method of steepest descent is an asymptotic approximation method for 
certain types of exponential integrals of the form 

It= [ g{z)e'f^'^dz (16) 
Jc 

where C is a contour in the complex z-plane and g{z) and f{z) are complex- 
valued analytic functions. The parameter t is taken to be real and positive, 
and we are interested in the asymptotic behavior of (jl6p as t — > oo with t > 

0. Laplace's and stationary phase methods arc just instances of this general 
procedure. The integral is dominated by the highest stationary points of /, 

1. e. if f{z) = u{x,y) + iv(x,y) with z = x + iy wc expect the integral to be 
dominated by points where u is maximum and v is constant. The only possible 
extrema for / are the saddle points, where f'{z) = 0. Since / is analytic, u and 
V satisfy the Cauchy-Riemann equation 

dx'^ dy"^ 

and from the maximum principle [53] we have that if > then < or 
vice versa. If zq is the saddle point, then we can deform the contour to C (by 
Cauchy's theorem) so that it passes through zq. From the Taylor expansion of 
f{z) about Zq we have 

f{z)-f{zo) + \.r{zo){z-zo)', 

where ~ means "is close up to additive error to". Then g{z) ^ g{zo), because 
for large t the main contribution to the integral comes from /. Then It becomes 

Jc 

Setting 

z-zo-re'^ and /"(zo) = |/"(zo)| e*^ 

it can be seen that 

It ^ 5(^o)e*^("«) / exp(it |/"(zo)| e''>+^^'^r^)e''^dr. 

Note that (j) is the angle of inclination of the oriented tangent to C at point 
zo, i.e. <f) = arg(zo) on C [23j. Choosing + 20 = tt, i.e., cj) = (tt — 9)/2 then 

/t - 5(zo)e*/(^")e"^ / e-5*k"(^°)l'-'dr 
Jc 
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and solving this Gaussian integrao yields 

/■^9(.o)«'"»>e'*(^)"%0r'). (17) 

The deformation of the contour chosen to make the integration Gaussian cor- 
responds to the steepest descent path from the saddle point, hence the name of 
the method [53] . Taking this path is not essential, other methods like stationary 
point and Perron's method take another path with similar results |23| . 



3.2 Asymptotic Approximation of the Walk on the Line 
3.2.1 Left Amplitude 

First the integral- form corresponding to equation (|14|) is solved. First, put the 
integral in the form of equation (jl6p by setting n ~ jt = n/t) and writing 

«r(70 = ;^ r J29Ak)e'f^^'^ (18) 

where 

f,{k) = log X J {k)-ik J, (19) 

(k) = Nj {k)(j (k) (^-m sin k ± Vb^ - sin^ fc) . (20) 

The saddle points 9j of fj{k) are defined by the equation 

sun \ ■ asin6'j 
fji^j) = -«7T . = = 0. 

/&2 -a2 sij^2^^. 



This equation has a solution at 



= ± arcsin ( — ] . (21) 



Also note that |Aj(0j)| — 1. Moreover 



fA^j) = -^iSj + log I I (22) 



^2 



and 

±(72 - 1)^5272 + a2(i_ ^2) 



I^iOi) = ^ i '^^'^ ' ' ■ (23) 



*Thc Gaussian integral or probability integral is given by e ^ dx = y/H. 
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Another solution to the equation f'{Oj) = is at — tt — 0j in the interval 
[— 7r,7r]. However, since f"{Oj) and /"(— tt — Oj) have similar behavior, the 
computations do not change. 

The contour is the real line in [— 7r,7r] and has no imaginary part, therefore 
(j) = aig Oj = in equation (fTT]) . 

Now using p7l) . the asymptotic expansion can be obtained 



27r 



1/2 



. 2^1-7^ 
27r|6| 



-ijOj t 



1/2 



i|72 - 11^5272 + a2(l_ ^2) 



+ ori). 



3.2.2 Right Amplitude 



Next is the solution of equation Following the same steps as above, write 
the integral as 



1 

2^ 



.tfAk) 



(24) 



where /, is defined in the same way as in (jl9p . and 



hj{k) = Nj{k)^j{k)bi 



Ak 



(25) 



Reusing the previous calculations for (equations ([21]) . ((22|) and ([SB]) ), the 
asymptotic expansion is 



a, 



1/2 



+ 0(t-i) 



±& + ^02(1-72)4-5272 

27r|6| 



1/2 



i|72 - 11^62^2 +a2(l_ ^2)1 



0(r 
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4 Closed-form Formulas and Convergence 
4.1 Formulas 

Approximate closed-forms for the amplitudes of the state of the walk on the line 
were given. Now the main contribution of this paper can be stated formally. 

Theorem 3. Let ^ n/t and a = e"^^ + e"^^ , h = e^'^^i - e"^^ . // the state 
of the walk is 



1*0 with \Mn)) 



then, 



±b{l - 7) 



where the terms Aj , Nj , and 9j are given by 



A, = 



±6+ v/a^(l-72)+62^2 



2yr 



27r|6| 



1/2 



i|72 - 11^6272 + 02(1-72) 



N, = 



—ia sin 6*1 ± 1/52 — gjjj2 ^ 



-1/2 



= a^(0) i^-iasinOj ± ^/62 - a2sin2 6', 

+ ao'(0)6*e-*''^ 
&7 \ 



sinfl., = ± 



with a^(0) and a^(0) as the initial amplitudes of the walk for n = 0, and 
'^cr(^) = CKo*("-) = for n 0. 

In a seminal work, Konno |121ll3j gave explicit expressions for the amplitudes 
of a [/ (2) coin using a discrete path integral method. However, these expressions 
were not in closed-form, as it is claimed in this work. 

In order to assess the quality of the approximation, figures [2] and [3] show a 
comparison between the probability distributions given by Theorem[3l and a nu- 
merical simulation of walks that start with an equal superposition of directions 
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0.10 
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-100 



-50 



50 



100 



Figure 2: Comparison between the probability distributions of numerical simu- 
lation (dark) and Theorem [3] (dashed) with ti = 1/2 and T2 = 0, t = 100, and 
initial state in equal superposition of directions. 



for different values of the parameters. It can be seen that the approximation 
gives some errors, but the asymptotic agrees with the simulation. The figures 
show that Theorem [3] is close to the real values of the probability distribution, 
in particular in the middle part of the plots. 

The errors in the approximation made by Theorem [3] can be computed from 
two parts, the Euler-Maclaurin formula and the steepest descent method [53]. 
Denote these errors by e and e respectively. Let Bi = J2l=o (D^i-r be a 
Bernoulli number [S], and let d e {<— , — >}. Then, the error for af{'^t) is J^j ^j-d+ 
Ej^d, where 



and 




where pj is either equation (|20p if d , or (|25p if d >. It can be seen that 
if we take m terms from each summation, ej^d = 0(2"™) and ej^d — 0{t^"^). 
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Figure 3: Comparison between the probability distributions of numerical sim- 
ulation (dark) and Theorem [3] (dashed) with ti = 3/4 and T2 = 1/2, t = 100, 
and initial state in equal superposition of directions. 



4.2 Convergence and Properties 

For quantum walks on the line and n-dimensional grids there exists weak conver- 
gence theorems [8]. In this section, we state the weak convergence of quantum 
walks on the line with phase parameters using these previous results. Then 
we show some applications of the convergence to compute the support of the 
probability density function. 

Theorem 4. Let Q = [— 7r,7r] x {1,2} be a probability space with probability 
measure /j. = \ {ipQ{k)\Xj {k))\'^ dk /2Tr for k G [— 7r,7r] and j = 1,2. Define a map 
/i : — > M such that for {k, j) € ft 

h{k,j)^h,{k)^{-iy ''""^ 



sin^ k + tan^ |(ti - T2) 

Let Xt be a position of the quantum walk at time t with distribution given by 
(QP, and Z be a random variable of with distribution ^. Then we have as 
t — >■ 00 

where =^ denotes weak convergenc^. 



''A sequence of random variables {Xi : i > 1} converges weakly to a random variable Z if 
lim„_>oo = Z, given that lim„^oo E[h{Xn)] = E[h{Z)] for all bomid continuous functions 
/i : M -)■ K. 
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Proof. Consider the theorem that states the weak convergence of quantum walks 
on the hue [H theorem 1]. Let Xj{k) be as in Lemma [2] Then 



A'(fc) 



Dividing this by Xj{k) we obtain 



-a sin k cos k sin k 



2\J\P- - a2 sin^ k 



-iX'jik) 



ai sin k 

— a? sin^(fc) 



Then, after some algebra and observing that | = e"/^ tan ■|(ti — T2), the the- 
orem follows. □ 

As an application of Theorem HI we can calculate the position of the two 
peaks of the walk for large time. 

Corollary 5. The limit distribution of Xf/t is concentrated on the interval 

H H" 
2 ' 2 

Proof. Theorem 3] have its maximum and minimum values for k = ±tt/2 and 
the corollary follows. □ 

The maximum probability of Pt{n) is found at the top of these two peaks, 
i.e., where n = ±|a|/2 [8]. Considering \n/t\ as the speed of the peaks, it can be 
seen that by setting ti = T2 it gets its maximum value, i.e. the fastest spreading 
of the walk. This corresponds exactly to an identity operator, and the walk does 
not mix at all inside the range of corollary [51 In order to get high speed and 
maximum randomness (i.e. the best mixing for positions inside the range) for 
Pt{n), we can set any value such that |ti — = 1/2. This implies that the 
support of h is in [-l/v^,l/v^]. In this case, the operator simulates exactly 
the probability distribution of a Hadamard operator [8] . 

As another application of ThcorcmUl we can compute the density function of 
the random variable Y = Xt/t in the asymptotic limit when t ^ 00. Following 
the steps of [8] for the Hadamard coin, we differentiate the quantity 



P{Y < y) 



E 



ke[-i^ ,iT\:hj(k)<y 



Mk) X,{k) 



2 dk 
2^' 



(28) 



which yields the density function 

fiy) 



\b\/2 



7r{y^ - l)v/(|a|/2)2-y2 



(29) 

for y G (— |a|/2, |a|/2), under the assumption of Im{a^ ■ a^^*) sin(Ti — r2)7r — 
and jag" I = Ic^iTI = which agrees with |T2l, ITS]. 
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5 Conclusions 



This paper presented a study of discrete-time quantum walks on the hne. A 
symmetric SU{2) coin operation was proposed and analyzed as a step towards 
an understanding of quantum walks. Using Fourier analysis and asymptotic 
approximation methods, we computed a closed-form formula for the amplitudes 
of the state of the walk. With this formula, we have a direct way to compute the 
amplitudes at any time step without recurring to time-consuming simulations 
or numerical integration. This also give us a complete characterization of the 
induced probability distribution of general quantum walks on the line. 

One important question that remains unanswered is the relation between 
Theorems [3] and HI Theorem [3] is based on the computation of saddle points of 
the high oscillatory kernel of Fourier coefficients. On the other hand. Theorem 
[4] is based on the method of moments (see [8] for details). A relation between 
these two density functions could set a common ground for the analysis of coined 
quantum walks. 
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